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Abstract 

A numerical program is presented which facihtates the computation of the full set of 
one-gluon loop diagrams (including ghost loop contributions) , with M attached external 
gluon lines in all possible ways. The feasibility of such a task rests on a suitably defined 
master formula, which is expressed in terms of a set of Grassmann and a set of Feynman 
parameters, the number of which increases with M. An important component of the 
numerical program is an algorithm for computing multi-Grassmann variable integrals. 
The cases M=2, 3, 4, which are the only ones having divergent terms, are fully worked 
out. A complete agreement with known, analytic results pertaining to the divergent 
terms is attained. 
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PROGRAM SUMMARY 

Title of program: DILOG 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ire- 
land 

Computer for which the program is designed on and others on which it has been tested: 
Personal Computer 

Operating systems or monitors under which the program has been tested: Windows 98 

Programming language used: FORTRAN 77 

Memory required to execute with typical data: 73 728 words 

No. of bits in a word: 64 

No. of processors used: one 

Has the code been vectorized or parallelized?: no 

No of bytes in distributed program, including test data, etc: 589 824 

Keywords: Computation of one-gluon loop Feynman diagram in QCD 

Nature of Physical problem 
The computation of loop diagrams in QCD with many external gluon lines is a time con- 
suming task, practically beyond reasonable reach of analytic procedures. In this paper we 
apply recently proposed master formulas for the computation of one-loop gluon diagrams in 
QCD with ah arbitrary number, M, external gluon lines. 

Method of solution 

The Grassmann variables in the master expressions and their properties are represented in 
suitably defined integer matrices, in order to carry out the Grassmnan integration first. The 
parametric functions are handled by an integer representation. The output of the program 
is the complete analytic result (tested for validity up to M = 4). 

Restrictions on the complexity of the problem: M must one exceed 4, though extension to 
higher values is simple. 

Typical running time: 25 seconds for M = 4 
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LONG WRITE-UP 
1. Introduction 

Quantum Chromodynamics (QCD) enjoys, in our days, universal acceptance as the fun- 
damental theory for the strong interaction. As a quantum field theoretical system, QCD has 
been extensively applied to situations in which its perturbative content provides a depend- 
able computational tool. It is, in fact, within the framework of this perturbative content that 
QCD has successfully confronted the quantitative description of the multitude of scattering 
processes, which probe strong interaction dynamics at high energies. Admittedly, the study 
of the non-perturbative domain of the theory offers intriguing challenges. Nevertheless, the 
immediate need to confront recent measurements coming from the HERA and Tevatron par- 
ticle accelerators as well as the expected ones, in the near future, from the LHC accelerator 
continues to put perturbative QCD (pQCD) to the forefront of theoretical activity^. 

Given the non-abclian structure of QCD, the (by far) most demanding component of the 
theory in relation to perturbative calculations is its gluonic, as opposed to its quark, sector^. 
In particular, perturbative computations involving Feynman diagrams with gluon/ghost 
loops become, to say the least, quite monstrous. During the last decade or so various meth- 
ods, aiming to expedite Feynman diagram computations in QCD, have been proposed whose 
basic feature is that they rely in a first, rather than the usual second, quantization approach 
to the formulation of the theory. Corresponding attempts have employed either strings [1-3] , 
or world-line paths [4-10] as their underlying basic agents. Within the framework of the 
latter case, two of the present authors, Refs [9,10], were involved in work which led to the 
formulation of a set of master expressions, that condense the multitude of all Feynman di- 
agrams entering a given configuration determined by the number of loops and the number 
of external propagators attached on them. To be more precise, the derived expressions go 
up to two loop configurations, nevertheless the 'logic' of the construction can be extended 
to loops of higher order. Even so, the analytical confrontation of a two loop situation with 
four 'external' gluon lines constitutes a challenging enough problem [11]. 

"'^Non-pcrturbativc input in the relevant computations enters in the form of 'initial' phenomenological 
information, with pQCD taking charge from thereon. 

^For that matter, this is more so the case for the non-perturbative domain of the theory. 



2 



The basic feature of the master expressions arrived at in Refs [9,10] is that they are 
furnished in terms of a set of Grassman and a set of Feynman variables. Once integrations 
over these two sets of variables arc performed one obtains the full result, i.e. the one which, 
for the given configuration, contains the contribution of all Feynman diagrams at once. It is 
obvious, even before laying an eye on these master formulas, that in order to put them into 
practical use, one should employ suitable computational methods for confronting them. It 
is the aim of this paper to present such a program, which will be applied to the one gluon 
loop case for two, three and four external gluonic lines. Given that one part of our program 
deals with the confrontation of multi-Grassmann variable integrals, which are entangled with 
expressions involving additional variables (also in line to be integrated over), it is hoped that 
it could find applicability to other situations, where Grassmann variables also make their 
entrance. 

Our paper is organized as follows. In the following section we present the battery of 
formulas, which are associated with the master expression corresponding to one gluon/ghost 
loop with M external gluon attachments in all possible ways. We intend to consider the 
cases M — 2,3 and 4, which exhibit divergent terms^, in addition to finite ones. In section 3 
we describe the structure of the program, while section 4 presents our results, with the first 
two cases being explicitly displayed. Finally, our concluding remarks are made in section 6. 

2. The one loop master formula 

Consider a configuration consisting of one gluon/ghost loop onto which M external gluon 
lines, with corresponding momenta pi ■ ■ -pu are attached (see Figure). According to Ref [9], 
the master expression, which summarizes the total contribution from all Feynman diagrams 
pertaining to this configuration is given by 

rr(Pi,---,PM) = --^/^(27r)V^) Trc(tr---^S^)7^ I dTT^-'x 



\n=l 



X 



■ 1 .1 

n L dUr. 
.n=M-'^ 



9{um, ■ ■ .,Ui)F^^^\ui, . . . ,M^;T)exp 



T Yl Pn-PmG{u 



-\-permutations , (1) 



^As expected, the aforementioned master expressions implicate the absence of divergent terms for M > 4, 
cf. Ref [9]. 
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where g is the couphng constant of the theory, the I'q , i — 1, . . . ,M are the SU {^) color group 
generators (in the adjoint representation) with Trc the trace over the color group, the Ui 
are Feynman parameters, the function 9 is specified by 
9{um, . . . , -Ui) = 0{um - um-i) ■ ■ ■ 0{u2 - ui) and 



F(^\m,...,UM;T) 



.n=M 



n d^nd^n 



X 



X exp 



(2) 



In the above equation ^'s are Grassmann variables, the are polarization vectors for 
the external gluons, $'^1 is the so-called spin factor entering the world-line description of 
QCD (see below), with Tr^ denoting trace with respect to Lorenz generator representation 
indices and the G(ii„, Um) are free propagators for the particle modes entering the worldline 
path integral description of QCD, in the context of its first quantized version (see Ref 
obeying the equation (s) 



-dndmG{Un,Um) = dlG{Un,Um) = G{Un,Um) = 2[6{Un,Um) - 1] , 



(3) 



with boundary condition 



dnG{Un, Um) = G{Un, Um) = sign{Un " Um) " 2{Un " Um) = -G{Um, Un) , (4) 

The explicit expression for the spin factor in terms of the set of parameters entering our 
expressions is (the J^i, are the Lorentz generators, in the vector representation) 



$W =Pexp 



M 



n] 



n=l 



HI/ 



^ M ^ M n2-l 

= Sf^u + lz{Jpa)^iw X] (t>pa{n) + {lzf{Jp2a2)ii\{Jpia^)\u 1] 0P2(T2 (^2)0piai (m) + ■ ■ ■ , (5) 

n=l n2=l ni=l 

where 



4 , 



n) = 2Un{elPn,u - e'lpn,,) + ^in+iCn+iUn{el-^'e: - 6:^'el)S{ur,+i - 



^p,v\ 



T 



Ur. 



(6) 



A point of note is the following: In the above expressions a specific time ordering has 
been chosen according to which index n + 1 comes immediately after index n, with ^m+i = 

^M+l = 0. 



3. Presenting the program structure 

The form of eq. (1) is considered as a sum of different elements. Every element is 
described with a line in three matrices: CI, C2 and C3. 

Matrix CI contains the Grassmann part of the form. It is a matrix with integer elements 
with dimensions {NCI x NC2)^. 

Matrix C2 has dimensions {NCI x 1) and contains real numbers. 

Matrix C3 represents the functions that accompany the Grassmann variables and its 
dimensions are {NCI x A^C3). The representation of the functions will be discussed in the 
"Functions Handling" subsection. 

3a. Grassmann Coding 

A Grassmann variable is represented in matrix CI by an integer number. This integer 
is equal to the index of the Grassmann number for the ^ type variables and equal to the 
opposite of the index of the Grassmann number for the ^ type variables. So 

^ n , In ^ -n . (7) 

A form that contains products of Grassmann variables is represented by a line in matrix 
CI. The elements in this line are put in the same order by which the Grassmann variables 
enter a given product. This line is accompanied by a line in matrix C2 which contains a 
real number and a line matrix C3 which represents the functions. For example the following 
coding would take place 

C2 C3 CI 

a ■ funci ■ func2 ■ 6666 ^ ( a )( /i /2 )( 1 -1 5 -5 ) , (8) 

where the integers /i,/2 represent the functions funci,func2 and will be discussed in the 
following subsection. 

The addition of two forms containing Grassmann variables (represented as ®) is carried 
as follows 

ai ■ funci ■ func2 ■ 66 + 012 ■ funcs ■ 6666 ^ 

C2 C3 CI C2 C3 CI 

/2)(2 -2 )e( q;2 )( /3 )( 3 -3 1 -1 ) ^ 
^The first number is the number of lines and the second the number of columns. 
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/ 



(9) 



C2 C3 CI 

«i V /i /2 U 2 -2 
V «2 J V /3 J 1^ 3 -3 1 -1 
The result of the addition of two forms represented with two sets of matrices C2, C3 and CI 
will be a new set of these matrices with new dimensions. If the dimensions of the first set of 

matrices are represented by A^'s and the these of the second set by M's, then the resulting 
dimensions of the final set will be 



C2 C3 CI C2 C3 CI 

(iVi X l)(iVi X N3){Ni X N2) © (Ml X l)(Mi x M3)(Mi x M2) ^ 

C2 C3 Cl 

[(iVi + Ml) X l][(iVi + Ml) X max(iV3, M3)][(iVi + Mi) x max(iV2, M2)] . 



(10) 



So the dimensions of the matrices which represent our forms is not constant and it is impor- 
tant for every part of the program to be known. Where it is necessary, the lines of matrices 
C3 and Cl are filled with zeros. The addition of two forms is carried out in the program by 
the subroutine "ADD" . 

The multiplication of two forms containing Grassmann variables (represented as <S)) is 
carried as follows 

(ai ■ funci ■ func2 ■ 66 + "2 • func^ ■ 6666) • • func4 ■ 66 + P2 ■ func^ ■ 66) ^ 




/l /2 /s 

h u 
\h h J 



(11) 



The result of the multiphcation of two forms represented with two sets of matrices C2, 
C3 and Cl will be a new set of these matrices with new dimensions. If the dimensions of 

the first set of matrices arc represented by A^'s and those of the second set by M's, then the 
resulting dimensions of the final set will be 

C2 C3 Cl C2 C3 Cl 

{Ni X l)(iVi X N3){Ni X N2) ® (Ml X l)(Mi x M3)(Mi x M2) ^ 
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C2 C3 Cl 

[(TVi • Ml) X l][{Ni ■ Ml) X (iVa + M3)][(Ari ■ Mi) x {N^ + M2)] . (12) 

The multiplication of two forms is carried out in the program by the subroutine "MULTI- 
PLY". 

Certain simplifications and rearrangements are carried out in the Grassmann and Func- 
tion Number matrices upon the completion of certain routines. This is done by subroutine 
"REARRANGE" which performs the following acts: 

• It checks if in the same line two equal Grassmann variables exist and if that happens 
it puts zeros everywhere in the hues of matrices C2, C3 and Cl. This accounts for the fact 

that en=en=o. 

• If a zero exists in one line of matrix C2, then the corresponding lines in all matrices 
are omitted, all the following lines are moved in front by one line and the line dimension is 
reduced by one. 

• It collects the Grassmann variables to the left of every line of matrix Cl, ignoring the 
zeros in between. If this reduces the dimension of the columns of matrix Cl, it performs the 
reduction. The same act is performed on function matrix C3. 

• It places the integers which represent the Grassmann variables in every line of matrix 
Cl in ascending order, according to their absolute value. Among two variables with the 
same absolute value the positive is placed to the left of the negative one. The sign of the 
corresponding value of matrix C2 is changed according to the changes made to the order of 
the Grassmann variables. 

• It places the number which represents the functions in every line of matrix C3 in 
ascending order. The set of functions which are not allowed to change their relative order 
are placed to the left. Functions represented by numbers Snm with n > m are entered as 
3mn. This practice, as well as the two preceding ones, are needed for comparing the different 
lines of the matrices. 

• If two lines in the matrix Cl and the corresponding lines in the matrix C3 are found 
equal, then the accompanying factors of matrix C2 are added. The second line is omitted in 
all matrices and their line dimension is reduced by one. 

• The lines of the matrices are also rearranged. The lines with all elements of matrices 
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CI and C3 equal to zero are entered first. 

The matrices C3 and CI are two-dimensional, but their dimension is not known a priori. 
This can cause problems when the matrices are fed from one subroutine to another. In 
every subroutine the dimensions of these matrices are pre-defined so that they will always 
be larger than what is required. Also, it is not always possible to define the same dimensions 
for these matrices in all the subroutines. The FORTRAN language stores the elements of 
a two-dimensional matrix in neighbouring memory places so that the first column is stored 
first, then the second, etc. When these places are read from another subroutine they will not 
result to the same matrix in the new subroutine, unless the matrix has the same dimensions. 
To avoid this problem the subroutine "LINEUP" is used, which rearranges the elements of 
a matrix that is defined with dimensions (A^i x A^2) so that they are stored exactly the same 
way as they would be stored if the dimensions of the matrix were (A^{ xN'^. In other cases we 
deal with this problem by using one- dimensional arrays to store two-dimensional matrices. 
The elements of the matrix a^, with dimensions (A^^i x A^2), are stored into the elements h of 
the array, so that the first column of the matrix is stored into the first elements of the array, 
then the second, etc., according to the correspondence 

ttij 6(j_i).iVi+i . (13) 

To evaluate eq. (1) we first code the following form in matrices C2, C3 and CI 

^ Y inlnimlmS^ ' £'"9„9^C(xi„, Um)5{Un+l - Un) ■ (14) 

This is done by subroutine "KNKM". Then the exponent of the above form has to be 
evaluated. For this purpose the subroutine "EXPONENTIAL" has been constructed. In 
this routine it is checked whether a part that does not contain Grassmann variables exists. 
If this is the case, the exponent of this part will only multiply all the lines of the rest. 
Then the exponent of the Grassmann containing part is calculated. This exponent will only 
contain M terms, since there are only M different Grassmann pairs $,n^n- Every term is 
constructed by multiplying the Grassmann part of "KNKM" as many times as it is needed. 
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Every newly constructed term is then added to the pervious ones. The result is stored to 
the output of subroutine "EXPKNKM" . 

The structure iijl^J ■ <^{ji) is coded in subroutine "FMN". The index n + 1 is evaluated 
from the specific time order which has been chosen. The exponent of "FMN" is found by 
subroutine "PEXPONENTIAL" . The trace of b^^ (the first term in the sum of eq. (6)) is 
accounted for by putting the factor 4 to a hne of matrix C2. 

The complete result is determined in subroutine "MULTIPLYALL" . In the output of 
"PEXPONENTIAL" a line with zeros in the domain of matrices C3 and CI but with the 
factor —2 in C2 is added, to construct the form TrjJ^^^ — 2. Then the result is multiplied 
with the result of "EXPKNKM" . Now the complete Grassmann part of the calculation is 
stored to matrix CI. The Grassmann integration is simple, since it amounts to keeping only 
the terms where all the Grassmann variables are present. Another simplification is made by 
neglecting the terms that contain only one function to be traced^ and which has zero trace. 
Our calculation at the end of "MULTIPLYALL" is stored only in matrices C2 and C3 since 
the Grassmann variables have been integrated out. 

3h. Function Handling 

In order to keep the matrices that represent our forms as small as possible, we have first 
carried out the Grassmann integrations in order to remove completely the Grassmann vari- 
ables. In fact, this is the reason we have coded the functions that accompany the Grassmann 
variables in matrix C3. Every line of this matrix contains a product of functions which are 
represented by integer numbers according to the labelling described in Table 1. 

In order for the representation of Table 1 to be valid it must be M < 9. The asterisk 
has the meaning that the corresponding functions are considered non-commutative and the 
order that they have been inserted is not allowed to change. Their Lorentz indices u are 
also considered indefinite, that is they are not defined at this point by the number which 
represents the function. So when in a line of matrix C3 functions of the type 2n or 3n are 
found together it is assumed that the Lorentz trace has to be performed, according to the 
exact order the functions are found in the respective line. For example if the numbers 21, 
^Function of the type 2n or 3n of the next subsection. 
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22 and 31 are found in one line of C3 (placed in this order and without any other function 
of the type 2n or 3n in the same line), then it is assumed that 



21 22 31 ^ (21)^^ (22)^, (31),^ ^ 

{i/2){Jafs)^p2(£lpi^/S - £jpi,a)(i/2)(J^5)p^2(£^p2,5 - £5P2,j) 

{t/2){J,,U{l/T){ele\-elel) 



(15) 



In this way we postpone the evaluation of the trace until after the Grassmann integration. 

After the Grassmann integration we are able to perform the trace in the surviving terms. 
This is accomplished by subroutine "TRACE". The information contained in matrix C3 
is divided into two matrices C3 and C6. The new matrix C3 contains only the following 
functions according to the labeUing 

Matrix C6 contains in its first column the exponent of T followed by internal products 
of the vectors e's and p's, according to the coding of Table 1. 

In order to arrive at that output, the maximum number of traces to be calculated is 
evaluated first, so that the dimension of the matrix that will contain the result is determined. 
The saturation of indices p, a in Jp^ is performed instantly, since 



{il2){Jp^)^,2{eyn,a - eyn,p) = -2 (£>„,. - ^X.^) , 



(16) 
(17) 



Then the indefinite Lorentz indices become definite. We begin by giving the indices 1,2 
to the first function that enters the product to be traced, then the indices 2,3 to the second, 
etc, until the last one which takes as second index 1. The representation used is 



3ni/ , Pn,u 2nu . 



(18) 



The products arc then represented as lines in two temporary matrices and every term in a 
sum uses a separate line, for example 



£lp2,2 - elp2,i 



/ 


l\ 


/ 


311 


222 




\ 




V 


312 


221 


/ 



(19) 
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The addition and the multiphcation of forms as the above is carried out in the same way 
as was done when we dealt with the Grassmann variables, until all the terms to enter the 
product to be traced are exhausted. In order to form internal products of the different 
elements, the last digit of every element is compared to the last digit of all the elements to 
the right in the same line. When two last digits are found equal an internal product between 
the corresponding elements is formed and represented according to the prescription of Table 
3. 

We have, also, taken into account the rule • p„ = and we have made everywhere in 
matrix C6 (but not in C3) the replacement pi — —p2 — pa . . . —pM, according to momentum 
conservation. 

The result that corresponds to one line of matrix C3 is put in one line of matrix C6. 
This contains a sum of different terms, each of which is a product of internal products 
accompanied by a factor. When the line of matrix C6 is read, an integer less or equal to 
2000 (which represents a factor) signals the beginning of the product of internal products 
and a another one signals its ending. For example the following translation is assumed from 
a line of matrix C6 

-1 3431 1 3433 3223 3122 2424 -3 3423 3332 ^ 

• £^ + £^ • £^ £^ • P3 £^ • P2 P4 ■ P4 - 3 £^ ■ P3 ■ . (20) 

4. Results 

At the beginning of its run the program asks the user to insert the number of external 
gluons M. Then it calculates all the possible combinations of time order and lists them for 
the user to see. Every combination is accompanied by an integer number and the user is 
expected to insert one of these numbers in order to choose the desirable combination. 

Then the calculation is carried out. Several files arc used during this calculation. The 
results after the Grassmann integration (matrices C2 and C3 in coding of Table 1) are 
written in file "TEMP2". The final result (matrices C2, C3 in coding of Table 2 and C6 in 
coding of integers) are written in file "TEMP4" . The result in function like output is written 
to file "TEMPS". The integers NCI, NC3 and NC6 are connected to the dimensions of the 
output matrices. These are C2(A^C1 x 1), C3(iVCl x iVC3) and CQ{NCl x NCQ). The 
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result is read as a sum of different terms. Every term is formed by reading a line of matrix 
C2, then the corresponding line of matrix C3 and then the corresponding line of matrix 
C6. When a lot of non-zero elements exist in some lines of matrix C6, these lines are folded 
beneath in order to permit easy reading. Then the whole sum is supposed to be multiplied 
by the function and integrated by the multiple integral which appear at the begging of the 
output. 

The number of terms at the output (NCI) grows rapidly when M increases. For M = 2, 
three terms exist (shown to Table 4a). For M = 3, 23 terms exist. These terms arc shown 
in Table 4b. For M — 4 the existing terms are 233. These include both finite and divergent 
terms. 

Our results have been compared with the analytic calculations for the divergent terms 
for M = 2, M = 3 and M = 4 [9] and have been found in total agreement. The structure of 
the program and the methods used for the computation of the master formulas permit the 
extension of the calculations for value of M higher than 4. 

5. Concluding remarks 

In this paper a computational algorithm has been presented for the successful compu- 
tation of the complete set of one-gluon loop Feynman diagrams with two, three and four 
external gluon attachments, on the basis of the master formulas derived in Ref [9]. It is 
hoped that the particular feature of the constructed algorithm, namely the ability to ex- 
pedite integrations over a multivariable set, a subset of which is Grassmannian, could find 
wider applications to analogous situations that may arise in other physical problems wherein 
Grassmann variables make their entrance. Within the context of the present application, it 
would be of interest to apply the particular algorithm developed in this work to the two- 
gluon loop M = 4 case, the corresponding master expressions for which have been derived 
in Ref [10]. As a first attempt, one could restrict the relevant computation to the divergent 
term associated with the M = 2 configuration and verify the consistency with second order 
corrections to the running coupling constant in pQGD. 
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Tables 



Function 


Number representing 




Function 


{i/2){Jp„)^,2{eyn,a - eyn,p) 

\ 1 / \ /-^ / r \ (L/-*- '"5'-' (J 1r / 


2n (*) 


(z/2)(J,,)^,(l/T)(£«+^£^ - £ri£^)5K+i - 


3n (*) 


M M 




exp{r^ ^ bn •PmG(^^n,^im)]} 


100 


n=l m=n+\ 




-Pm dnG{Un,Um) 


2nm 


{l/2T)e^ • e^dr,dmG{ur,,u^) 


3nm=3mn 



Table 1. The representation of functions through integers before the Grassmann integration 

in matrix C3. 



Function Number representing 

Function 
- li") 3n 

M M 

exp{T^ J2 [Pn- PmG{Un,Um)]} 100 
n=l m=n+l 

■PmdnG{un,Um) 2nm 
• e"^dndmG{um Um) 3nm=3mn 

Table 2. The representation of functions through integers in matrix C3 at the output. 



Function 


Number representing 




Function 


Pn ' Pm 


2n2m 


-Pm 


3n2m 




3n3m 



Table 3. The representation of functions through integers in matrix C6 at the output. 
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M= 2 

NC1,NC3,NC6= 3 3 4 
TIME ORDER= 1 2 

-(pi**2/2)*g**(2)*dl4(pl+p2 )*TrC(tGa2 tGal )* 

inf 1 u2 

I dT I du2 1 dul th(u2,ul ) 


4.0 EXP(all) T**(-l) +2 e2.el p2.p2 

2.0 EXP(all) el.p2dlG(ul,u2) e2.pld2G(u2,ul) T**(-l) 

2.0 EXP(all) el.e2dld2G(ul,u2) T**(-2) 



Table 4-a. The function like output of the program for M = 2 and for the time order 
ui < U2. The Memoradum is similar to Table 4b- 
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M= 3 

NC1,NC3,NC6= 23 4 41 
TIME ORDER= 12 3 

-(pi**2/2)*g**(3)*dl4(pl+p2+p3 )*TrC(tGa3 tGa2 tGal )* 

inf 1 u3 u2 

I dT I du3 I du2 | dul th(u3,u2,ul ) 




;.0 EXP(all) 



T**( 0) 



4.0 


EXP (all) 


cl.p2dlG(ul,u2) 








0) 


4.0 


EXP (all) 


el.p3dlG(ul,u3) 








0) 


4.0 


EXP (all) 


e2.pld2G(u2,ul) 








0) 


4.0 


EXP (all) 


e2.p3d2G(u2,u3) 








0) 


4.0 


EXP (all) 


e3.pld3G(u3,ul) 






fp**l 


0) 


4.0 


EXP (all) 


e3.p2d3G(u3,u2) 








0) 


2.0 


EXP (all) 


el.p2dlG(ul,u2) 


e2.pld2G(u2,ul) 


e3.pld3G(u3,ul) 




0) 


2.0 


EXP (all) 


el.p2dlG(ul,u2) 


c2.pld2G(u2,ul) 


c3.p2d3G(u3,u2) 




0) 


2.0 


EXP (all) 


el.p2dlG(ul,u2) 


e2.p3d2G(u2,u3) 


e3.pld3G(u3,ul) 


pi**! 


0) 


2.0 


EXP (all) 


el.p2dlG(ul,u2) 


e2.p3d2G(u2,u3) 


e3.p2d3G(u3,u2) 


fp**l 


0) 


2.0 


EXP (all) 


el.p3dlG(ul,u3) 


e2.pld2G(u2,ul) 


c3.pld3G(u3,ul) 


pi**! 


0) 


2.0 


EXP (all) 


cl.p3dlG(ul,u3) 


e2.pld2G(u2,ul) 


c3.p2d3G(u3,u2) 


pi** 1 


0) 


2.0 


EXP (all) 


cl.p3dlG(ul,u3) 


e2.p3d2G(u2,u3) 


c3.pld3G(u3,ul) 


p>**l 


0) 


2.0 


EXP (all) 


el.p3dlG(ul,u3) 


e2.p3d2G(u2,u3) 


e3.p2d3G(u3,u2) 


pi**^ 


0) 


8.0 


EXP (all) 


dl(u3-u2) 






pi**l 


-1) 


8.0 


EXP (all) 


dl(u2-ul) 






pi**l 


-1) 


2.0 


EXP (all) 


el.p2dlG(ul,u2) 


e2.e3d2d3G(u2,u3) 




pi**l 


-1) 


2.0 


EXP (all) 


el.p3dlG(ul,u3) 


e2.e3d2d3G(u2,u3) 




pi**l 


-1) 


2.0 


EXP (all) 


e2.pld2G(u2,ul) 


cl.c3dld3G(ul,u3) 




p>** 1 


-1) 


2.0 


EXP (all) 


e2.p3d2G(u2,u3) 


el.c3dld3G(ul,u3) 




p>** 1 


-1) 


2.0 


EXP (all) 


e3.pld3G(u3,ul) 


el.e2dld2G(ul,u2) 




p>** 1 


-1) 


2.0 


EXP (all) 


e3.p2d3G(u3,u2) 


el.e2dld2G(ul,u2) 




pi**! 


-1) 



-1 


e3.p2 


e2.p3 


el.p2 






+1 


e3.e2 


el.p2 


p3.p2 






+1 


c3.c2 


cl.p2 


p3.p3 






-1 


e3.p2 


e2.el 


p3.p3 






+1 


eS.el 


e2.p3 


p2.p2 






+1 


e3.el 


e2.p3 


p2.p3 






-1 


e3.e2 


el.p3 


p2.p2 






-1 


c3.c2 


cl.p3 


p2.p3 






-1 


e3.el 


e2.p3 


p3.p2 






+1 


e3.p2 


e2.p3 


el.p3 






+2 


c3.p2 


c2.p3 


-2 


c3.c2 


p3.p2 


+2 


e3.p2 


e2.p3 


-2 


e3.e2 


p3.p2 


-2 


e3.p2 


el.p3 


+2 


e3.el 


p3.p2 


+2 


e3.el 


p3.p3 








-2 


e3.p2 


el.p3 


+2 


e3.el 


p3.p2 


+2 


e3.el 


p3.p3 








-2 


e2.p3 


el.p2 


+2 


e2.el 


p2.p2 


+2 


e2.el 


p2.p3 








-2 


c2.p3 


cl.p2 


+2 


e2.el 


p2.p2 


+2 


e2.el 


p2.p3 









-2 

+2 



e3.p2 
eS.el 



e2.el 
e2.p3 



+2 e3.el 
-2 e3.e2 



e2.p3 
el.p3 



Memorandum: 
pi**2 

dl4(pl+p2+p3 ) 

TrC(tGa3 tGa2 tGal ) 
inf 



5'*^ (Pl +P2 +P3) 

Trc{tlHlHl^) 



th(u3,u2,ul ) 
EXP (all) 

d2d3G(u2,u3) 
dl(u3-u2) 



S(w3,«2,ui) 

M M 

exp{r^ ^ [p„ 

n=\. m.=n-\-l 

d2dzG{u2^uz) 

(5(W3 - U2) 



Pm G(u„,Um)]} 



el.p2 



dlG(ul,u2) 9iG(ui,U2) 
T**(-l) T-i 



Table 4b- The function like output of the program for M — 3 and for the time order 

Ui < U2 < Us- 
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Figure Caption 

Figure Illustration, for M = 4, of the classes of one-gluon loop Feynman diagram (right 
side of arrow) which are simultaneously accommodated by the corresponding master formula 
depicted on left side of arrow. 
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